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Abstract. We put forward a method of constructing discrete coherent states for n qubits. 
After establishing appropriate displacement operators, the coherent states appear as displaced 
versions of a fiducial vector that is fixed by imposing a number of natural symmetry 
requirements on its Q function. Using these coherent states we establish a partial order in the 
discrete phase space, which allows us to picture some n-qubit states as apparent distributions. 
We also analyze correlations in terms of sums of squared Q functions. 



1. Introduction 

Coherent states (CS), first introduced by Schrodinger (T), are of paramount significance for 
modern physical theories, as they are quantum states that follow classical trajectories. In 
quantum optics, CS were popularized by Glauber J2][3) for the description of correlation 
properties of a single-mode radiation field, where the Weyl-Heisenberg group emerges as a 
hallmark of noncommutativity |4j . 

Although generalizations along several directions have been considered (see references 
(5]-[7] for reviews), nowadays it seems indisputable that the pioneer work of Perelomov JSJ 
paved the way to extend the notion of CS to quantum systems with dynamical symmetry 
group G. In this approach, CS appear as orbits of a certain fiducial state under a unitary 
irreducible representation of G acting in the corresponding Hilbert space. 

This fiducial state is chosen to have maximal isotropy subgroup H, which ultimately 
leads to "maximal classicality". Under these circumstances, the displacement operators, 
transforming CS among themselves, are labeled by points in the manifold M = G/H. For 
many physical models, M can be equipped with an irreducible symplectic structure |5}|T2j, so 
it can be considered as the phase space of a classical dynamical system. In other words, there 
is a one-to-one correspondence between CS and points of the classical phase space. 

The crucial point of this construction is that the Hilbert space is irreducible under 
the action of G. This is especially clear for the symmetric representations of the unitary 
groups SU («), when the classical phase space is a (2« — 2) -dimensional sphere [13]. For the 
discrete counterparts, even in the physically (but not mathematically!) simple n-qubit case, 
the symmetric subspace is only a very small portion of the whole 2" -dimensional Hilbert 
space. Nonetheless, one can define a natural set of CS, constructed in a similar way as in the 
continuous case: acting on a fiducial state with discrete displacements; i.e., unitary operators 
labeled by elements of two discrete sets p4) . These two sets can be organized in a discrete 
2" x 2" grid, on which a specific discrete geometry (including symplectic operations) can be 
introduced, so that such a grid turns out to be a bona fide discrete phase space 1 15-23). 



Although the points of the discrete phase space label again CS, there is still an essential 
difference with the continuous case: the choice of the fiducial state. For continuous symmetry 
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groups, the standard choice corresponds to an extreme state of the representation, such as 
the vacuum or the lowest/highest weight state. In the discrete case, the nature of the unitary 
displacements prevents such a simple notion and different possibilities have been discussed 
sofarj24l|25). 

In this paper, we take the fiducial state as a spin CS and impose that its associated 
Q function fulfills reasonable symmetry conditions. This not only solves the problem, but 
allows us to use the system of CS to impose a partial order in phase space, which helps 
to recognize states pictured as distributions. Finally, we briefly speculate about detecting 
quantum correlations through the sum of squared Q functions. 



2. Discrete phase space 

A qubit is a two-dimensional quantum system, with Hilbert space isomorphic to C 2 . It is 
customary to choose two normalized orthogonal states, {|0), |1)}, as a computational basis. 
The unitary operators 

a z = |0)<0|-|l)(l|, a.- = |0)(l| + |l)(0|. (2.1) 

generate the Pauli group 3?\ under matrix multiplication J26| . 

For n qubits, the Hilbert space is the tensor product C 2 ® • • • ® C 2 = C 2 .A compact way 
of labeling both states and elements of the corresponding Pauli group & n consists in using 
the finite field F2« (27). This can be considered as a linear space spanned by an abstract basis 
{9i, . . . ,9,,}, so that given a field element a (henceforth, field elements will be denoted by 
Greek letters) the expansion 

n 

a = £a i i , a, eZ 2 , (2.2) 
i=i 

allows us the identification a (a\, . . . ,a n ). Moreover, the basis can be chosen to be 
orthonormal with respect to the trace operation (the self-dual basis); that is, 

tx(e i e j ) = 8 ij , (2.3) 

where tr(of) = cc + cc 2 + ... + a 2 " ' , which actually maps ¥2" H> Z2. In this way, we associate 
each qubit with a particular element of the self-dual basis: qubit, 0, . 
The generalized Pauli group is generated now by the operators 

Z a = £>(aA)|A)(A|, X p =Z\X + p)(X\. (2.4) 

X X 

Here the additive characters % are defined as x( a ) = exp[/7Ttr(a)] and |A) is an orthonormal 
basis in the Hilbert space of the system. Operationally, the elements of the basis can be 
labeled by powers of a primitive element (i.e., a root of the primitive polynomial), and read 
{|0), \a),..., Iff 2 " -1 = 1)}. One can verify that 

Z a Xp= X (aP)XpZ a , (2.5) 

which is the discrete counterpart of the Weyl-Heisenberg algebra for continuous variables |4|. 



The operators ( 2.4 1 can be factorized into tensor products of powers of single-particle 
Pauli operators. This factorization can be carried out by mapping each element of F2» onto an 
ordered set of natural numbers according to 

Z a = of 1 ® . . . ® of" , Xa = a b x l ® . . . ® o h » , (2.6) 
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where a,- = tr(a0,) and /?, = tr(j30,) are the corresponding expansion coefficients for a and j3 
in the self-dual basis. Moreover, they are related through the finite Fourier transform |28| 



^ = 4=I*(AA')|A><A'|, (2.7) 
v2" a ^ 



so that X M = PZp J?. 



We next recall |20] that the grid defining the phase space for n qubits can be appropriately 
labeled by the discrete points (a, j3), which are precisely the indices of the operators Z a 
and Xp: a is the "horizontal" axis and j3 the "vertical" one. In this grid one can introduce 
the concept of straight lines (also called rays), which possess the same properties as in the 
continuous case. It is worth noting that the monomials labeled by points of the same ray, 
{(a,fia)} commute with each other, so that one can establish a correspondence between 
eigenstates of such commuting sets and states (actually bases) in the Hilbert space p0[ . 

Following our program, we introduce the set of displacements 

D(a,P) = ^ a ^Z a X p> (2.8) 

where <I>(a,j3) is a phase required to avoid plugging extra factors when acting with D. One 
can immediately check that 

D(a,p)D*(a,P) = l, D f (a,p)=D(a,P), (2.9) 

so they are unitary and Hermitian. They also constitute a complete trace-orthonormal set 

TT[D{a,P)D{a',p)]=2 n 8 aa ,8 w . (2.10) 

These operators act multiplicatively on the monomials \2A\ , thus shifting phase-space points 
according to 

(a,/3) D( ^' ) (a + a',j3+j3'), (2.11) 
which justifies their designation. 

Finally, for later purposes, we touch on a pair of symplectic operations (z- and x- 
rotations) that transform rays into rays according to 

P^Pl - ZaX^a , Q v X p Q\ - Z v pX p . (2.12) 

The symbol indicates equality except for a phase. Both P^ and Q v are unitary operators, 
with [PjjjXy] = [Q v ,Zn] — 0, and can be written as 

^=£^|A)<A|, e v = £c^ v |A><A|, (2.13) 

x x 

where |A) are the eigenstates of Z a and \X) of Xp. The coefficients c% v fulfill the recurrence 
relation 

CX+a.y =c a , v CA,vZ(vaA), C 0)V = 1 , (2.14) 

whose explicit solution can be found in reference [29 1 and it is unimportant for the rest of the 
paper. 

We associate an eigenstate \yo) of Z a with the horizontal axis and immediately obtain 
that the state associated with the ray j3 = jJ.CC is Pn\yo), while the vertical axis is associated 
with J?\Yo) |30") . Any other straight line, parallel to a given ray, but crossing the axis j3 at 
the point j3 = q, corresponds to the state X^P^^Iq). 

Using phase-space coordinates these z- and jc-rotations can be interpreted as 

(a,j3)^(a,j3 + M«), («,J3) H (a + vj3,j3). (2.15) 

It is clear that these two transformations are conjugate each other. 
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3. Discrete coherent states for n qubits 

According to the conventional approach, we define discrete CS | Of, J3 ), labeled by phase-space 
points (o£,/3), as the displacements of the fiducial state l^o) fl9) : 

|a,j3)=D(a,j3)|*F ). (3.1) 

The state I^Po) can be chosen in several ways p4) . Here, for reasons that will be apparent 
soon, we take | as a product of identical qubit states: 

|¥o) = |#,<P>i®---® (3.2) 

where 

\V, ( p) J =e^hm(^) Mj + e-Wco&fl) \ 0)j> (3 . 3) 



and the angles (#,<p) parametrize the Bloch sphere. The state (3.2 1 is invariant under 
permutation of the qubit indices and thus can be expanded as (24) 



with £, = e"? tan#/2. The basis {\k,n) : k = 0, . . . ,n} are the Dicke states 



M = V t PtWh- ■ -WkWm ■ ■ ■ |0)„) ■ (3.5) 



Here, {^} denotes the complete set of all the possible qubit permutations. 
In field notation, the state ( |3.4[ > can be compactly expressed as 

where the function Ii(k) counts the number of nonzero coefficients kj in the expansion of K 
in the field basis (see the appendix for a brief account of its properties). 

Finally, notice that one might think in imposing that the states |^) are eigenstates of the 
Fourier operator J3l][32) (much as the vacum is for continuous variables). Since J^ 2 = 1, this 
is tantamount to 

^) = ±IS), (3.7) 

which leads to two possible candidates E,± = ±\/2 — 1 and all the qubits pointing in the same 
direction. However, we prefer to follow an alternative route to fix the possible values of £, . 

3.1. P -function 

Let us look for an expansion of the density matrix of the form 

P = £p(a,j3)|a,j3)(a,j3|, (3.8) 

which is the analogous to the Glauber-Sudarshan P-function for continuous variables |33][34j. 
It is not difficult to check that the function P(a, j3) may be recast as 

P(a,j8)=Tr[pA(a,j8)], (3.9) 
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where Tr (with capital T to distinguish it from tr, which is the trace in the field) stands for the 
ordinary trace operation in Hilbert space. The kernel A(a,/3) reads 

A («>£) = ^LxW + fiy) ^\D(r,8\^)- l D(Y,d). (3.10) 

This clearly shows that the P-function is nonsingular only when (E, |D(y, 5)|<^) _1 exists. 
Since the state 14) is factorized into single-qubit states, we have 

(^l^^l^ocfj^tDlafcrfl^ 1 )), (3.11) 

where gi,di are the expansion coefficients of y and 5 in the self-dual basis and 

= -f=^=(|0) + §|1)). (3.12) 



Using equation ( 3.6 1 one can find 

(I;\d(y,$M) ' 



1_|£|2\ t*(r)-*(5)+*(r+«)]/2 



\ [A(«)-A(r)+*(r+«)]/2 / e_s* \ [ft(r)+fi(5)-/>(y+5)]/2 



,1 + I^IV U + IS 

This obviously rules out some values of 4 for which P is singular. 

As an illustrative example, consider the case of a single qubit. Then, one finds that 



(3.13) 



P{a,b) 4 



1 , l + IS 



2 



(3.14) 



where now a,/? € Z2. This function is singular when 4 is real, imaginary, and when |4| 2 = 1, 
i.e. on the equator of the Bloch sphere. In particular, the eigenstates of the discrete Fourier 
transform, when 4± — V2± 1, do not lead to the faithful expansion on CS for qubits, contrary 
to what one could expect. 



3.2. Q-function 

In our search for determining the values of 4, we next look at the Q-function, defined in 
complete analogy with its continuous counterpart, namely 

Q p (a,P) = (a,p\p\a,P), (3.15) 

which satisfies 

£<2 p (a,j3)=2". (3.16) 

a.p 

Let us impose the maximal symmetry conditions admissible on the Q-function for the fiducial 
state 

CI .- <2|£) (of,/J) is symmetric under axis permutations. 

C2 .- The values of <2|£)(ot,j3) can be obtained from Q^(a,0) by z- andx- rotations. 
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Since Q^(a,fi) = |(^|D(a,j3)|^)| 2 , and using equation (3.13 », one can easily find out 
that the condition CI [Q^(a,fi) = Q^(j5, a)] imposes the following restriction on £, (this 
is one of the possible restrictions, but the other ones are just symmetric reflections): 



t, = ^-\/l +cos 2 $ — cos exp(7#) , 
so that the 2-function takes the form 

cos # 



-n/2 < & < tc/2, 



(3.17) 



%)(«,j3): 



\f \ +cos 2 ■& 



2/j(ce+j3) 



sin# 



V 1 +cos 2 ■& 



[h(a)+h(P)-h(a+P)] 



.(3. 



To fulfill the condition C2, we first require that (a, 0) contains all the possible values 
of (2|£)(05,/J). To this end, we note that 

2*W / cos ^ x2A(a) 



C2|{)(k,k) 



sin# 



Vl +cos 2 # 



C2|{)(a,0) 



Vl + cos 2 # 



(3.19) 



so that the only possibility (apart from symmetric reflections) is that sin# = cos#, i.e. 

(3.20) 



# = tc/4. Consequently, equation (3.18 i takes the simple form 

j \ \h(a)+h(P)+h(a+P)] 



which explicitly fulfills Q\^(a, a) = Q\%)((X, 0). 

Using the properties of the function h (see appendix A), one can infer that for any ordered 
pair (0f,j3) there is always a field element K, given by 



K= Y,(ai + bi-aibi)0i, 



(3.21) 



{9j} being the self-dual basis, such that (ct,j3) = (k",0). This means that each value 
of (a, j3) can be obtained by rotating Q||)(K',0) according to the following protocol: 

(k,0) 4 (jc,|UJc) ^ (k+k:juv,k:ju), (3.22) 
where the rotation parameters are fi = j5 K~ 1 and v = (a + K") /3 1 . 



To conclude, we observe that, for a single qubit, the P-function (3.14i for the value £ in 

a+b Tr(po y ) 



equation (3.17 1 (at # = 7C/4), becomes 

i Vi r 



P(a,l 



■I) 



(3. 



4 4 L (-l) fo Tr( P (T z ) + (-irTr(p(T. 

which is the most uniform possible. Note, in passing, that this maximum uniformity could be 
also adopted as a reasonable criterion to fix the value of S, . The associated unit vector clearly 
reflects this uniformity 

n = ma x \^)^\a y ^),^\a l \^)) = 1/73(1,1,1). (3.24) 



4. Ordering points in the discrete phase space 

The very simple form of the Q-function for the fiducial state \£, ) in the previous section allows 
to introduce a partial order in the n-qubit phase space. Indeed, since any Qi^\(a,P) can be 
obtained by rotations from (fc,0) = 3~ h ( K \ we can order the points on the horizontal axis 
according to the values of the /i-function: < h(a) = k < n. The C\ = k ^"!_ k y_ elements 
which correspond to the same value of the /z-function remain disordered inside the strip with 
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Figure 1. Ordered Q-function for the fiducial state £ ) . 

the fixed value of h(a) = k. This automatically arranges the rest of the phase-space points 
according to the symmetry property and the construction ( |3.22| >. 

In figure 1 we plot the Q-function for the fiducial state |4 ) for 5 qubits using this ordering. 
Explicitly, the order of axis is chosen as: 0, (J 6 , a 26 , a, a 7 , a 2 , a 16 , a 27 , a 8 , a 17 , a 18 , a 3 , 
a 19 , 20 , e\ 6> 28 , e 10 , a 21 , a 22 , a 29 , a 23 , a 24 , a\ a 11 , a 5 , a 12 , a 25 , a 13 , a 14 , a 30 , a 15 , 
a 31 . The irreducible polynomial used is x 5 +x 2 + 1 =0, and the self-dual basis chosen for 
F 2 s is 0! = a 3 ,0i =ff 5 ,0i =o n ,0 1 ^a 22 ,e x ^a 2A . 

The shape of the Q-function presents a hump localized at the origin. Due to the 
covariance under displacements, Q\yS)( a iP) — Q\%)( a + Y>P + £>), the ordering should be 
applied to the pairs (a + 7, j3 + 5), but not to (a , j3 ) itself. In this sense, one cannot properly 
say that Q\^g) ( a >P) has a hump located at (7, 8) if we keep the previously established order. 
Nonetheless, it is clear that due to the functional form of 2|y,5) ( a > P) the elements of the field 
can be easily rearranged (using the summation table) in a such way that the corresponding 
hump becomes centered at (7, 8) and has a symmetric form. In figure 2 we plot the Q-function 
for the - qubit CS |0 10 , 10 ) = D(9 U \ 1O )|^) according to such a prescription. 

It is also worth observing that the Q-function of an arbitrary state I 1 ?) can be written 
down as 

<2 w (a,j8)=£P(a + 7,j3 + 5)<2| 5) (7,5), (4.1) 

i.e. as an smearing of the P-function. In particular, the order established by the /i-function 
helps to visualize the superpositions of several discrete CS as spatially separated humps in 
phase space. In figure 3 we plot the Q-function for a superposition of two CS for 5 qubits, 
\V) oc +D(9 3 \ 31 )|£)), and the ordering is the same as for 2|g) • We can clearly observe 
two humps with a residual symmetry. 
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Figure 4. Ordered symmetrized Q-function of the fiducial coherent state | £ ) displaced to the 
center of the phase-space for eight qubits. 



As a final remark we may note that the distribution still can be re-ordered in a more 
symmetric form just distributing the points with the same value of h(a) on both sides of the 
principal peak. Although the number of such points is not always even, for a large number of 
qubits the distribution corresponding to the CS is practically symmetric, as it can be seen from 
figure 4, where we plot the 2-function for the fiducial state corresponding to 8 qubits. 

in the limit 



It is worth comparing the form of the 2-function for the n-qubit CS (3.1 



n ^> 1 with the CS resulting from taking as the fiducial state an eigenstate of the discrete 
Fourier transform pT) . In the latter case the (^-function tends to a Gaussian shape fT9j , 
while in our approach it has a step form modulated by a decreasing function f(k), which 
along the axes Z (j3 = 0), X (a = 0) and Y (a = j3) has an exponential form: 



/(*) = S*o + 3-* £ 

m=0 



m+1 



e U-£c? -e It- £c; 

V r=0 J \ r=0 j 



(4.2) 



where 9 (k) is the Heaviside step function. 



5. Detecting correlations in n-qubit systems 

By construction, the discrete CS are factorized states, so the qubits therein do not exhibit 
correlations. For symmetric states, the correlations are frequently measured using the concept 
of spin squeezing [35-40], comparing the fluctuations of some definite operator with the 
standard quantum limit, given by the spin CS. Nevertheless, for nonsymmetric states similar 
criteria do not work well |41 1 for the choice of the measured operators becomes nontrivial. 
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Figure 5. Q-function for the state XORi^l^). 



To study correlations in nonsymmetric n-qubit states we apply a criteria proposed in 
reference [42] to quantify polarization fluctuations. According to this approach, we compute 
the sum of squares of the Q-fxmc\\on: quantum correlations make such a sum lesser than the 
corresponding one for a CS. In fact, for the fiducial state ||) we have 

/ , n lh(a)+h(P)+h(a+P)] „ l / -.\[2a i +2b i -2a i b i ] /A \n 

l^mheQ) =n e (3) =m • (5.D 

a,p a,p \ 3 s i=l<i,A=0^ ' \ ' 

To check the method, we consider a simple way to induce correlations between qubits: 
the application if XOR M gates, where the pair (p,q) indicates the qubits on which the operator 
is applied, namely 

xoRp !q \ai,...,a p ,...,a q ,...,a n ) = \ai,...,a p ,...,a 9 +a p ,...,a n ) . (5.2) 

For a correlated state I 1 ?) = XOR p q D(ji, v)|4) the sum of Q 2 curiously does not depend on 
the form of the displacement D(jj.,v) and gives 



which is smaller than ( 5.1 1. In the same vein, the application of k XOR gates between different 
particles (i.e., now |¥) = XORp i)ft ...XORp 1)9l £>(/x,v)|€), with Pl ^ qi ^ ... ^ p k ^ q k ) 
keeps decreasing the sum: 

5*><«-»-(£)*Gr- 
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Figure 6. Q-function of the squeezed state S a i \£,). 

Similarly, the application of sequences of XOR gates to the fiducial state also leads to 
decreasing values of the Y*Q 2 { a T&)- This effect can be clearly seen in figure 5, where the 
Q-function for the state XORi 2^) is plotted. One can observe that the heights of the Q(a,j5) 
are smaller, so that the distribution initially localized at the origin sparse over a substantial 
part of the phase space. 

To induce correlation between all the qubits in a regular way one can apply the squeezing 
operator |21 25] 

S C = J»<£A|, (5.5) 

which acts on the Z a and Xp as the scaling transformation 

S^XpS^. =Xpr , SrZ a Sl =Z„j-i , (5.6) 



much as in the continuous case. The action of ( |5.5| l on a CS can be formally expressed as 

^ > = (TT^?^ V) ' (5 - 7) 

and implies that the initially factorized state (|3.2|i is transformed into 



where 



d .i = L fijmCm ■> c m = tr(Cftn) , fijm = tx(QiOjO m ) . (5.9) 

i,m= 1 
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The squeezing operator correlates all the qubits in a generic CS (3.1 1 and the degree of 
such correlation depends both on % and £. For example, in the 5-qubit case the operator 
S a 7 correlates qubits in the initial CS with | = e m l A (V? — l) /s/2 in the most efficient way 



according to the criteria (5.1 1, and its action on the fiducial state |^) is 

1 

x %{ki+k2+k3+k 4 }^{k l+ k 2 +k 5 }\ k ^ k ^ k ^ k ^ k5 ^ (5 1Q) 

and {• + •} means sum mod 2. In figure 6 we plot the 2-function of the state S a i \£,), where it 
can be observed that the initial distribution is spread out over practically all the phase space. 

6. Conclusions 

We have developed a method for constructing discrete CS from the symmetry conditions for 
the Q-function of the fiducial state. This has allowed us to order the points in the discrete 
phase space. Besides, we have applied a criterion for the detection of quantum correlations to 
the discrete case and have shown that it can be useful for «-qubit systems. 
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Appendix A. Some properties of the function h(a) 

The function h(a) is defined as the number of nonzero components in the expansion of a field 
element a in the self-dual basis {0, }, that is 

n 

where a,- = tr(a0,). Note that < h(a) < n. The basic properties we need in this paper are 
the following: 

£ X(ao t )= £(-ir=n-2fc(a), 

i=i i=i 



(A.2) 



h(a+P) =h(a)+h(P)-2Y,aibi 

i=i 



where b\ = tr(/30, ). The second of these equations follows from the equality 

h(a + P) = {a\ + bx} + {a 2 + b 2 } + . . . + {a„ + b n } . (A.3) 

Here {• + •} denotes again the sum mod 2 and verifies 

{a,- + bi} = a t + bi - laibi . (A.4) 
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